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GENERATION OF THE INVARIANT COEFFICIENTS 
OF THE CHARACTERISTIC POLYNOMIAL 
FOR AN n x n MATRIX 


1. Introduction 

Often in the theories of numerical stability, one is seeking roots to a char- 
acteristic polynomial, which in the case of the predictor with iterative correction 
method of numerical integration, are eigenvalues of a matrix whose elements 
depend on the coefficients used in the integration process. 

It is the purpose of this document to explicitly display the characteristic 
polynomial in terms of the elements of the characteristic matrix and code the 
same for future numerical analysis. 


2. Analysis 

Let a tj be the elements of an n x n square matrix; and let the characteristic 
polynomial in X. be defined as the eigenvalue determinant 

D(\) = | 3ij -XS..| 

The explicit representation of this polynomial is 

n 

»w = Lh^ 

4=0 


where the coefficients, b^, can be obtained via Taylor’s expression as 


h - 1 d ^ D 


\= o 


The plan for obtaining an explicit representation of b^ will be to differentiate 
D (a.) from its explicit representation as a sum of products of the elements 


[a 


k, j 


- \ K J 


l 



Consider the representation 


*-L 

P { j • | i = 1 , 2 * * n} 

n J 1 1 ’ 


n 



k = 1 


where the sum is over all permutations of the n elements jj > j 2 • • • 3 n *> e j * s 
±1 depending on whether the permutation is even or odd, and i k is the ordered 
sequence {1, . . . n}(i k = k); and the set { j. | i = 1, 2, . . . n} are the same set of 
numbers, {1, 2 . . . n}, and j k is the k th element of the permutation. 

The coefficient b Q is obtained by setting A. = 0 in D(A ). 

n 

b n = T e. | | [a. . ] 

0 3 I | l k* 3 k 

P n { J i I i = 1 » 2 * * n} k = 1 


In accordance with the plan, we differentiate D (A.) with respect to A. and 
obtain 


dD 

dA 



J2 e i s k\j k , TT [a i k i k 

P n {j i |i=l,2--n} Vk* 


k 8 . 


Now only those permutations for which j k » = k' will be non-zero, and since j k » 
does not explicitly appear in the product, we may write 

~=- T" e. ~TT [a. . -AS. . ] 

d A / I / 1 3 1 k ,3 k *k * 3 k 

k -1 P n _i { J i I i = l. 2* • n^k'} k = 1 

where no\y the set j. excludes k’ and the sequence i k is the ordered permutation 
of the set of numbers { j . } . Meaning to e . is obtained by recognizing that the 
ordered permutation has e . = +1. 
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We proceed now to take other derivatives, each time we obtain a negative 
sum over new k’s different from previous k’s and the set j excludes the numbers 
taken by the set of k’s in the sums. Also, the product sequence i k becomes the 
ordered sequence of the j k set. We obtain for example the second derivative as 


d 2 D 
d\ 2 


- (- 1) 5 


E IE 




P n _ 2 { j A | i = l * •*n^k^ 1 \ k^ 2) ) 


n- 2 


k = 1 


[a. . -KS. . ]• 


H J 


k J k 


^'Jk 


In general, the ^ th derivative becomes 


^ = (-!)< 
dk* 


£ £ - ,£«., 
k (1) = i k (2) Vk (1) k ( ^Vk .•••k (1) 

n -'t 

Y e.TT [a. -\S. ]. 

, i M l k ,J k H’Jk 

p n _ / £ {j j | i= 1 . * * ‘ n^k (1 ) * * * k^ } k = 1 


Now, since only the set {k (1) . . . k^ } }, and not its ordered sequence enters 
into the summand, many terms in the multiple sum on the k's will give rise to the 
same contribution. In particular, since {k^ 1 ) . . . k^} are distinct, there are 
exactly V. ways to obtain the sequence k (1 ) . . . k ( 5 from a particular combi- 
nation of i numbers. As a consequence, the multiple sum on k’s can be written 
as V. times a sum over all combinations of l numbers taken out of the n numbers, 
1, 2 . . . n. Since the permutation P n _^ are over the remaining elements, we 
could just as well choose all combinations of h - ^ elements out of the n numbers 
and permute over this set. Explicitly, we obtain 


d^D 

dk^ 


(- 1)^ l\ 


L 

^nU 




k 8 , 


where 





is the sum over all combinations of n -i numbers out of the set of 
numbers 1, 2 ... n; 
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y * is the sum over all permutations of the n-Z numbers obtained from 


P n-^ 


e . 


each of the preceeding combinations. 

is +1 if the sequence in the permutation is even (The ordered sequence 
has e. = 1) and is -1 if the sequence is odd. 

And the set i k is the ordered sequence { j k } of the n -l numbers from . 

The explicit representation of b^ is obtained by setting k = 0 and we obtain 




(- 1 ) 


l 


L Z> 

c lu p n u 


n -'t 


k = 1 


[ a i 3 J 

k ' J k 


3. Computer Program 

The program implementing this formulation for b^ consists of a MAIN 
routine and two subroutines, COMBIN and PERMUT. MAIN reads in a single 
integer N (2 - N - 9), the order of the matrix; then, on the basis of N creates 
the coefficients b a , b 2 , . . ., b N . (b 0 is not created since this is merely the 
determinant of the matrix [a . .] ). The coefficients are created and printed in 
literal, FORTRAN format. A symbolic FORTRAN program similar to the 
printed version can be had as a JOB CONTROL LANGUAGE option. 

For an orderly generation of the coefficients, MAIN relies on the two sub- 
routines COMBINE and PERMUT. COMBIN produces, on (^) successive calls 
all combinations of the n items 1, 2, . . ., n taken m at a time. The combinations 
are generated in ascending order, (i.e., 1, 2, 3, < 1, 2, 4 < 1, 2, 5 < 1, 3, 4 etc.) 

PERMUT produces on m! successive calls all permutations of m items and 
the sign of the permutation. PERMUT requires that the items to be permuted be 
in ascending order on the first call. With this requirement we can define each suc- 
cessive permutation as positive or negative obviating the time consuming chore 
of checking the result of a permutation. 

As a matter of interest and as an indication of the practical limits on using an 
explicit representation of these coefficients, Table I gives the number of terms 
in each coefficient for various orders of n. 
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Appendix A 
Sample Output 


LITERAL FORM OF 
THE INVARIANT COEFFICIENTS OF 
THE CHARACTERISTIC POLYNOMIAL FOR A 

4X4 

MATRI X 


F ( 4) = +1 


F( 3 ) = 
F ( 2 ) = 

* 

$ 

F( 1) = 

* 

* 

* 

* 

* 


* 

* 

* 

* 


+A ( 1 » 1 ) + A( 2,2 )+ A(3» 3)+A{ 4,4) 

+A(3,3)*A(4,4)-A( 3, 4 )* A{ 4, 3) +A ( 2 , 2 ) * A ( 4, 4 ) - A ( 2 ,4 > *A (4,2) 
+A(2,2)*A(3,3)-A(2,3)*A(3,2)+A( 1,1)*A(4,4)-A< 1,4 )*A(4,1) 
+A ( 1,1 >*A< 3,3 ) -A ( 1,3)*A{ 3, 1)+A< 1, 1 )*A<2,2 ) - A ( 1,2 )*A(2,1) 
+A ( 2, 2)* A (3, 3 )*A<4, 4)-A(2,2 )*A( 3,4 )*A<4,3) 

+A ( 2, 4) * A( 3 » 2 ) *A (4 , 3 ) - A( 2 , 4) * A ( 3, 3) * A ( 4, 2 ) 

+A(2, 3)*A(3,4)*A(4,2)-A(2»3)#A( 3 , 2 ) *A ( 4, 4) 

+A< 1,1 )*A(3,3 )*A(4 ,4)-A< 1, 1)*A(3,4)*A<4,3) 
+A(1,4)*A(3,1)*A(4,3)-A( 1,4)*A( 3,3)*A<4,1 ) 

+A < 1,3) *A< 3,4)*A<4 , 1 )-A< 1, 3)*A(3, 1 )*A(4,4) 

+A (1,1)*A( 2,2 ) *A (4, 4)-A ( 1, 1 )* A( 2,4)*A(4,2) 

+A( 1 , 4 ) * A ( 2 » 1 ) *A (4 , 2 ) - A ( 1,4)*A(2,2)#A(4, 1 ) 

+A < 1, 2)*A( 2,4 )*A<4, 1 )-A< 1,2 >*A( 2,1 >*A<4, 4) 

+A( 1,1 )*A< 2, 2)*A< 3, 3)-A( 1, 1)*A(2, 3)*A( 3,2 ) 

+A(1,3)*A( 2,1)*A(3,2)-A< 1,3)*A( 2,2)*A(3,1 ) 

+A( 1,2)*A<2,3 )*A(3, 1 )-A( 1,2)*A(2, 1)*A(3,3) 
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LITERAL FORM OF 


THE INVARIANT COEFFICIENTS OF 
THE CHARACTERISTIC POLYNOMIAL FOR A 


5X5 


MATRIX 


F( 5)* 

+1 



F (4) - 

+A( 1 

1 )+A( 2 

, 2 ) +A ( 3 

F( 3 ) = 

+A (4 

4)*A< 5 

,5 )-A (4 

* 

+A ( 3 

3)*A<4 

,4 ) -A ( 3 

* 

+A (2 

2)*A<4 

,4)-A(2, 

* 

+A( 1 

1 )*A< 5 

, 5 ) -A ( 1 

* 

+A ( 1 

1)*A<3 

,3)-A(l, 

F (2) = 

+A ( 3 

3 ) * A ( 4 

,4 ) *A ( 5 

* 

+A (3 

5)*A(4 

,3 ) *A ( 5 

* 

+A (3 

4)*A<4 

, 5 ) *A ( 5 

* 

+A ( 2 

2)*A(4 

,4 ) *A ( 5 

* 

+A ( 2 

5 ) *A ( 4 

,2)*A(5 

* 

+A (2 

4)* A( 4 

♦ 5 ) *A ( 5 

* 

+A ( 2 

2)*A<3 

,3 ) *A ( 5 

* 

+A (2 

5)*A(3 

,2 )*A(5i 

* 

+A ( 2 

3)*A< 3 

» 5 )*A ( 5 

* 

+A ( 2 

2)*A( 3 

,3 )*A<4, 

* 

+A (2 

4)*A(3 

» 2 ) *A (4 

* 

+A (2 

3)*A( 3 

,4 ) * A ( 4 

* 

+A ( 1 

1 ) * A ( 4 

,4 ) *A ( 5 

# 

+A ( 1 

5 )* A( 4 

, 1 ) * A ( 5 1 

* 

+A( 1 

4) *A ( 4 

t 5)*A( 5 

* 

+A(1 

1)*A(3 

,3 ) *A ( 5 1 

* 

+A( 1 

5 ) *A ( 3 

, 1 ) *A ( 5 

* 

+A (1 

3 )* A< 3 

, 5 ) * A ( 5 1 

* 

+A( 1 

1 ) * A ( 3 

,3 ) *A (4 

* 

+A (1 

4)*A<3 

,1)*A<4, 

F(2)- 

F ( 2 ) +A ( 1 

3 ) *A ( 3 

,4 ) *A (4 

* 

+A (1 

1 )* A( 2 

, 2 ) * A ( 5 1 

* 

+A ( 1 

5>*A (2 

,1 )*A< 5 

* 

+A (1 

2 )* A( 2 

, 5 ) *A ( 5 i 

* 

+A< 1 

1 )*A< 2 

, 2)*A(4 

* 

+A ( 1 

4)*A< 2 

, 1 ) * A ( 4 

* 

+A ( 1 

2 ) *A ( 2 

,4)*A(4 

* 

+A ( 1 

1 )* A( 2 

, 2 ) *A ( 3 1 

* 

+A ( 1 

3>*A( 2 

, 1 ) *A ( 3 


+A ( 1 

2 )* A( 2 

,3 ) #A (3 i 


1.3) +A(4,4)+A(5,5> 

► ,5)*A(5,4)+A(3,3)*A<5,5)- 
),4>*A(4,3>+A(2,2>*A<5,5)- 
!,4)*A(4,2)+A( 2 , 2 ) *A ( 3, 3 )* 

,5)*A(5,1)+A(1,1)*A<4,4)- 
., 3)*A( 3, 1>+A< 1,1)*A(2»2)- 

1.5) -A(3,3)*A<4,5)*A<5,4) 
i, 4)-A( 3,5)*A( 4,4 )*A(5,3) 

1.3) -A(3,4>*A<4,3)*A(5,5> 

1.5) -A(2,2)*A( 4,5)*A<5,4) 
>»4)— A(2»5)*A(4»4)*A(5»2 ) 

1.2) — A(2,4)*A( 4 , 2 ) *A ( 5 , 5 ) 

1.5) -A(2,2)*A(3,5)*A<5,3> 

1.3) -A(2,5)*A( 3,3>*A(5,2> 
i»2)-A(2»3)*A(3»2)*A(5»5) 
>,4)-A(2,2)*A( 3 ,4 ) *A ( 4, 3 ) 

»3)-A(2»4)*A(3,3)*A(4»2) 
>» 2 )-A( 2, 3 )* A( 3 » 2 ) *A ( 4, 4 ) 

, 5 ) -AC 1, 1) #A ( 4, 5 ) * A ( 5, 4 ) 

, 4)-A( 1, 5)*A( 4,4 )*A(5,1 ) 

, I >— A ( 1,4) *A ( 4, 1 ) * A ( 5, 5 ) 

, 5 )-A ( 1, 1)*A( 3 ,5 ) *A ( 5 , 3 ) 
,3)-A(l,5)*A(3,3)*A(5, 1) 

, 1 ) — A ( 1,3)*A< 3,1 )*A(5,5) 
,4)-A(l,l)*A(3,4)*A(4, 3) 

► , 3)-A( 1,4)*A( 3,3)*A<4,1 ) 

, 1 ) — AC 1,3>*A(3,1)*A(4,4) 

i, 5 )— A< 1, 1)*A( 2 » 5 ) *A ( 5, 2 ) 

> , 2 ) - A( 1, 5>*A(2,2)*A(5, 1 ) 
i, 1 )-At 1,2 >* A( 2,1 )*A(5,5) 

, 4 ) -A( 1, 1) *A ( 2,4 ) *A ( 4, 2 ) 

► ,2)-A( 1,4)*A( 2,2)*A(4,1 ) 

, 1 )-A( 1 * 2 ) *A ( 2, 1.) * A ( 4, 4 ) 

l, 3)-A( 1, 1)*A( 2 » 3 ) *A ( 3, 2 ) 
(* 2 )-A( 1, 3) *A (2, 2 )*A ( 3, 1 ) 
i,l) —A ( 1 , 2 ) * A( 2,1 )*A(3,3) 


•A( 3 » 5 )*A( 5,3) 
•A ( 2 , 5 ) *A (5 ,2 ) 
A( 2 , 3 )*A< 3,2) 
• A ( 1,4) *A (4,1) 
•A ( 1,2 )*A( 2,1) 
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F ( 1 ) = +A ( 2 
+A ( 2 
+A ( 2 
-A ( 2 
-A ( 2 
-A (2 
+A ( 2 
+A ( 2 
+A<2 
-A ( 2 
-A ( 2 
-A ( 2 
+A ( 1 
+A ( 1 
+A ( 1 
-A ( 1 
-A ( 1 
-A ( 1 
+A ( 1 
+A ( 1 

F ( 1 ) = F ( 1 ) +A ( 1 
-A ( 1 
-A ( 1 
-A ( 1 
+A { 1 
+A ( 1 
+A ( 1 
-A ( 1 
-A ( 1 
-A ( 1 
+A ( 1 
+A ( 1 
+A ( 1 
-A ( 1 
-A ( 1 
-A ( 1 
+A ( 1 
+A ( 1 
+A ( 1 
-A ( 1 

F ( 1 ) = F ( 1 )- A ( 1 
-A ( 1 
+A ( 1 
+A ( 1 
+A ( 1 
-A ( 1 
-A ( 1 
-A ( 1 
+A ( 1 
+A ( 1 
+A ( 1 
-A ( 1 
-A ( 1 
-A ( 1 
+A ( 1 
+A ( 1 
+A ( 1 
-A ( 1 
-A ( 1 
-A ( 1 


, 2 ) - A ( 3 *3 )4A<4 , 
*2>4At3,5)*A<4, 
»2)*A( 3,4 )4A<4 , 
,5 )* A( 3,2 )*A (4, 
,5)4A(3,4)4A(4 , 
,5 >*A< 3,3 )4A<4, 

* 4 ) * A ( 3 * 5 ) 4A (4 i 
, 4)- A< 3,3 )4A<4, 
,4)*A(3,2)*A(4, 
, 3)4 At 3 ,4 ) 4A <4, 
,3)*A< 3»2>*A<4 , 
,3>*A(3,5)*A<4, 
,1 ) =:= A ( 3,3)4A<4 , 
, 1 )*A< 3,5 )*A(4, 
,1 ) * A ( 3,4>*A<4 , 
,5>*A(3,1)*A(4» 

* 5 ) *A ( 3 *4 ) 4A (4 , 
,5)*A(3,3)*A(4, 
, 4 ) - A ( 3,5>*A<4 , 
,4)4A(3,3)*A(4, 
,4)*A( 3, 1 >*A(4 , 
, 3)*A( 3,4)*A(4, 
, 3) *A( 3, 1)*A(4 , 
, 3)4A( 3,5 ) * A ( 4 , 
,1 )*A< 2,2>*A(4 , 
, 1)*A( 2, 5 ) * A ( 4 , 
,1 )*A< 2,4)*A<4 , 
,5 ) * At 2, 1 ) *A ( 4, 
, 5 ) * A ( 2 *4 ) 4A (4 , 
» 5 )4 At 2, 2 ) * A ( 4 , 
, 4 ) * A ( 2* 5 ) *A (4 , 
, 4 ) ^ A ( 2,2 ) *A (4, 
, 4) *A ( 2, 1 ) 4A (4 , 
, 2 ) * A( 2 , 4 ) 4 A ( 4 , 
,2)*A(2, 1>*A<4 , 

* 2 )4 At 2,5)4A<4, 
,1 >4A< 2,2)*A( 3, 
, 1 )4 At 2, 5 ) * A ( 3 , 
,1 )*A( 2,3 )*A( 3, 
» 5 )* At 2, 1 ) -A (3 , 
,5)*A( 2,3 )4A< 3, 
, 5 ) * At 2, 2 ) =*= A ( 3 , 
, 3 ) 4 A ( 2 , 5 ) 4A ( 3 , 
, 3)* At 2,2 ) -A (3, 
, 3)*A(2,1)*A(3, 
» 2 )* At 2,3 >*A(3, 
,2)4 At 2 , 1 ) 4 A ( 3, 
, 2 ) 4 At 2, 5 ) * A ( 3 , 
,1 ) *A < 2, 2) 4A ( 3, 
, 1 )4A( 2,4 )4A(3, 
,1 ) * A ( 2 , 3 ) *A ( 3, 
,4) 4 At 2, 1 ) 4 A ( 3 , 
,4)4A(2,3)*A(3, 
, 4)4A( 2,2 )4A(3, 
, 3 ) 4 A ( 2 ,4 ) 4A ( 3 , 
, 3 ) 4 A ( 2,2 ) 4 A ( 3 , 
, 3 ) 4 A ( 2 » 1 ) *A ( 3, 
, 2 ) 4 At 2,3 ) 4 A ( 3 , 
,2 )4A( 2, 1 ) 4A ( 3, 
, 2 ) * At 2,4 ) 4 A ( 3 , 


4 ) 4 At 5, 5) -A ( 2 
3 ) 4 A ( 5 , 4 ) — At 2 

5 ) 4 A ( 5,3) -At 2 
3 )*A ( 5, 4 ) + At 2 
2 )4 At 5, 3 ) +A ( 2 

4) 4A( 5,2 )+A( 2 

2 ) 4 At 5, 3) -A (2 
5 )4A ( 5, 2 ) - At 2 

3) 4A( 5, 5) -A ( 2 

5) 4A(5,2)+A( 2 

4) 4A(5,5)+A(2 

2 >4A< 5,4)+A( 2 

4 ) 4 At 5, 5) -A ( 1 
3)*A(5,4)-A( 1 

5 ) 4A( 5, 3) -At 1 

3) 4A(5,4) + A( 1 
1 )4A( 5, 3) +A ( 1 

4 ) 4 A ( 5 , 1 >+A( 1 

1 )4A( 5 , 3 ) - A ( 1 

5 ) 4A( 5, 1 ) - At 1 

3 ) 4 At 5, 5) -A ( 1 
5 ) 4 A < 5, 1 ) + At 1 
4 ) 4 At 5 , 5 ) +A ( 1 
1 )4A( 5,4)+A( 1 

4) 4fi| 5, 5) -At 1 

2 )4A( 5, 4)- At 1 

5 ) 4 At 5, 2 ) -A ( 1 
2 ) 4 A ( 5,4)+A( 1 

1 ) 4 A< 5, 2 ) +A ( 1 

4) 4 A t 5, 1 ) + At 1 
1 >4A( 5, 2 ) -A ( 1 

5 ) 4A( 5, 1 ) - At 1 
2 )4 At 5, 5) -A ( 1 
5)*A(5,1)+A( 1 

4 )4 At 5, 5) +A I 1 

1 ) 4 A < 5 , 4 ) + At 1 
3 ) 4 At 5, 5) -At 1 

2) 4A(5»3)-A( 1 
5 )4A( 5, 2) -At 1 

2 ) 4 A ( 5 , 3)+A( 1 

1 ) * A ( 5 , 2 ) + A ( 1 

3) 4A< 5, 1 ) +A( 1 

1 ) 4A< 5,2) -A ( 1 
5 )4A( 5, 1 )— At 1 

2) 4 At 5, 5 ) - A ( 1 

5 ) 4 At 5, 1 > + A( 1 

3 ) 4 A ( 5, 5) +A t 1 

1 ) 4A< 5, 3)+A( 1 

3 ) 4 At 4 , 4 ) - A ( 1 

2) 4A( 4, 3 ) -At 1 

4 ) 4 At 4, 2 ) -A ( 1 

2 )4 At 4, 3 ) + At 1 
1 ) 4 A ( 4, 2 ) + A ( 1 

3 ) 4 A ( 4, 1 ) + A( 1 

1 ) * At 4,2) -At 1 
4) 4 A t 4, 1 )-A( 1 

2) 4 At 4, 4) -At 1 
4 ) 4 A ( 4, 1 ) + At 1 
3 ) 4 A ( 4, 4 ) + A ( 1 
1 >4A( 4, 3 ) + At 1 


, 2 ) 4A t 3, 3 )4 A ( 4 
,2)4A( 3, 5 )4A( 4 
,2)4A(3,4)4A(4 
, 5 ) *A ( 3 , 2 ) 4 A ( 4 
,5)4A(3,4)4A(4 
,5 )4A( 3, 3)4A(4 
, 4 ) 4 A t 3, 5 ) 4 A ( 4 
,4 ) 4A( 3, 3 ) # A ( 4 
, 4 ) 4 A ( 3, 2 ) 4 A 1 4 
, 3 ) 4A ( 3, 4 ) 4A ( 4 
, 3 )4 A ( 3, 2 ) 4 A ( 4 
,3)4A(3,5)*A(4 
, 1 ) 4 A ( 3, 3 ) 4 A 1 4 
, 1 ) 4A( 3, 5)4A(4 
, 1)4A<3, 4)4 A (4 
, 5 ) 4A( 3, 1 >4A< 4 
, 5)4At 3, 4)4 a(4 
, 5 ) 4A ( 3, 3 ) 4 A ( 4 
, 4 ) 4 A ( 3, 5 ) 4 A ( 4 
,4 )4A( 3, 3 ) 4 A ( 4 
, 4 ) 4 A ( 3, 1 ) 4 A ( 4 
» 3 ) 4A ( 3, 4 ) 4 A ( 4 
, 3 ) 4 A ( 3, 1 )4A(4 
, 3 ) 4A ( 3, 5 ) 4 A ( 4 
, 1 )4A<2,2 ) 4 A ( 4 
, 1 ) 4A ( 2 , 5 ) 4A( 4 
, 1 )4A(2,4)4A<4 
, 5 ) *A ( 2 , 1 ) 4 A ( 4 
, 5)4A(2,4)4A<4 
,5)4A(2,2)4A(4 
,4)4A(2,5)4A(4 
,4)4A(2,2)4A(4 
, 4 ) 4 A ( 2 , 1 ) 4 A ( 4 
, 2 ) 4A ( 2 , 4 ) 4 A ( 4 
, 2 ) 4 A ( 2, 1 )4 A ( 4 
, 2 ) 4A ( 2 , 5 ) 4A ( 4 
,1)4A(2,2)*A(3 
, 1 )4A( 2, 5)4A( 3 
,1)4A(2,3)4A(3 
, 5 ) 4A ( 2 , 1 )4A(3 
,5)4A(2,3)4A(3 
,5)4A(2,2)4A(3 
,3)4A(2,5)4A(3 
, 3 ) 4A < 2 , 2 )4A(3 
, 3 ) 4 A ( 2, 1 ) 4 A ( 3 
,2)*A(2,3)4A(3 
, 2 ) 4 A ( 2 ♦ 1 ) 4 A 1 3 
,2>4A(2,5)4A(3 
,1)4A(2,2)4A(3 
, 1 )*A<2, 4)*A( 3 
,1)*A(2,3)4A(3 
,4 ) 4A ( 2 , 1 ) 4 A ( 3 
,4)*A(2,3)4A(3 
,4 ) 4A( 2,2 )4A( 3 
, 3 ) 4 A ( 2 , 4 ) 4 A ( 3 
, 3) 4A( 2,2 )4A( 3 
,3)4A(2,1)*A(3 
,2)4A(2,3)4At3 
, 2 ) 4 A ( 2 , 1 ) 4 A ( 3 
,2)*A(2,4)4A(3 


,5)4A (5,4) 
,4 )4At 5, 3) 
, 3 ) *A ( 5 ,5 ) 
,4 )4A( 5,3) 
, 3 ) *A ( 5 ,2 ) 
,2 )4A( 5,4 ) 
» 3 ) *A ( 5 » 2 ) 
,2 )4A( 5, 5) 
,5)4A (5,3) 
,2 )4A( 5,5) 
,5)4A (5,4) 
,4 )4A( 5, 2) 
,5)4A (5,4) 
,4 )4A( 5, 3) 
,3 )4A (5,5 ) 
,4 )4A( 5,3) 
, 3 ) *A ( 5 , 1 ) 
,1 )4A (5,4) 
,3 )*A(5, l) 
, 1 )4A( 5,5) 
,5 )*A (5,3) 
, 1 )4 At 5, 5 ) 
, 5 ) 4A ( 5 ,4 ) 
,4 )4A( 5,1) 
, 5 ) 4 A ( 5 ,4 ) 
,4 )4A( 5, 2) 
,2 )4A (5,5) 
,4 )4A( 5, 2) 
, 2 ) 4A ( 5 , 1 ) 
, 1 )4 At 5,4) 
, 2 ) 4A ( 5 , 1 ) 
, 1 >4A( 5,5) 
,5)4A (5,2) 
,1 )4A ( 5,5) 
, 5 ) *A ( 5 ,4 ) 
,4 )4At 5, 1 ) 
,5)4A (5,3) 
,3 )4A( 5,2) 
,2 )4A (5,5) 
,3 )4A< 5, 2) 
,2 )4A (5,1) 
, 1 )4A< 5,3) 
,2 )4A (5,1 ) 
, 1 )4 A ( 5,5) 
,5)4A (5,2) 
,1 )4A( 5 , 5 ) 
, 5 ) 4A ( 5 , 3 ) 
, 3 )4At 5,1) 
,4)4A (4,3) 
,3 )4A( 4,2) 
,2 >4A (4,4) 
, 3. )4 A t 4 , 2 ) 
,2 )4A (4,1) 
, 1 >4A( 4, 3) 
, 2 ) 4 A ( 4 , 1 ) 
, 1 )4A( 4,4 ) 
,4 )4A (4,2) 
, 1 )4 A ( 4 , 4 ) 
,4 )4A (4, 3) 
,3 )4A( 4, 1 ) 
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Appendix B 
Program Listing 

PURPOSE 

TO PRODUCE IN LITERAL FORM (FORTRAN FORMAT) THE INVARIANT 
COEFFICIENTS OF THE CH ARACT BU S T IC POLYNOMIALS FOR N-SOUARE 
MATRICES WHERE N MAY RANGE FROM- THREE TO NINE 

***** **************************************************** ************** 


N = DIMENSION OF N-SOUARE MATRIX 


N 1 = NUMBER OF CALLS TO COMBIN 
= BINOMIAL COEFFICIENT (N I) 

N2 = NUMBER OF CALLS TO PERMUT 
= I FACTORIAL 

N3 = NUMBER OF TERMS PER LINE (I FACOTRS PER TERM) 

= INTEGER 8/1 

N4 = NUMBER OF FACTORS PER LINE 
= I*N3 

N 6 = NUMBER OF ITEMS IN COMBINATION (N THINGS TAKEN Nf> AT A TIME) 
= N-I 


M 

= MAXIMUM 

1 VALUE 

OF I 

Ml 

= N-l 
COUNTS 

CALLS 

TO COMBIN 

M2 

COUNTS 

calls 

TO PERMUT 

M3 

COUNTS 

NUMBER 

OF TERMS 

M 5 

COUNTS 

NUMBER 

L I NE S ON 

I 

= NUMBER 

OF FAC 

TORS IN A 


COUNTS 

UP TO 

N— 1 = M 

K 

COUNTS 

UP TO 

I (NUMBER 


STORED FOR CURRENT LINE 
CURRENT PAGE 

TERM, NUMBER OF ITEMS BEING PERMUTED 
OF ITEMS IN PERMUTATION ARRAY) 


V COMBINATION ARRAY 

LA PERMUTATION ARRAY LA ( 1 , 1 ) TO LA(I,1) ARE THE RIGHT HAND 
INDICES FOR THE LITERAL OUTPUT FACTORS 
X INITIAL PERMUTATION I T BIS X(l) TO X ( I) ARE THE LEFT HAND 
INDICES FOR THE LITERAL OUTPUT FACTORS 
W THE PRINT ARRAY 


ITEMS 


NUMBER OF ITEMS BEING COMBINED AT A TIME 
N-I+l 


C#**** **************************************************** ************** 

C 

IMPLICIT INTEGER(A-Z) 

COMMON/NP/IS (9) ,LA( 9,9 ) , I TE MS, I FIRST 
DIMENSION SIGN(3),V(9),W(3,8),X(9) 

DATA SIGN/'-', '*','+•/ 
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100 READ! 5 * 100 0* EN 0=999) N 

I F( (N.LT.2>.0R.N.GT.9)G0 TO 999 
M = N — 1 

WRITE(6,20 00)N,N 
WRITE! 6,2010) 

WR I TE ( 7 , 2 600) N 
WRITE! 6,21 00 ) N 
WRITE (7,2100)N 

WRITE! 6,22 00) M, ( ( SIGN! 3) , I , I > , I = 1, N) 

WRITE(7 ,2200)M, ( ( SIGN! 3) , I , I) , 1=1, N) 

N 1 = N 
N2 = l 
M5=2 

DO 200 1=2, M 
N6=N-I 

Nl= ( Nl* ( N 6+ 1 ) ) /I 

N2=N2*I 

N3=8/I 

N4=I*N3 

ASSIGN 160 TO L 
V! 1 ) = 0 
M3=0 
M4= 0 

IF! M5+- ( N 1 *N2 ) /N 3 . G T . 60 )M 5 = 60 
DO 200 M 1 = 1 , N 1 
CALL COMB IN!N,N6,V) 

DO 110 J= 1 ,N 
110 LA! J,1)=J 

DO 120 J= 1 ,N 6 
120 LA! V (J ) , 1 ) = 0 
K = 0 

DO 130 J =1 , N 

IF<LA( J,1 ).EQ.O) GO TO 130 
K=K+ 1 

LA! K , 1 ) = LA ( J, 1 ) 

X! K) =LA ( J, 1) 

130 CONTINUE 
I FIR ST= 1 
ITEMS =1 

DO 190 M2 = 1 ,N2 
CALL PERMUT 
S=SIGN( I S < 1 )+2) 

DO 140 K= 1 » I 
W( 1 ,M3*I +K ) =S 
W(2,M3*I+K)=X(K.) 

W(3,M3*I+K)=LA(K ,1) 

140 S=SIGN( 2) 

M3 =M 3+ 1 

I F( <M3 .LT.N3) .AND. ( N2-M 2+N1 -Ml . N E .0 ) )G0 TO 190 
IF! M3.LT .N3)N4=I*M3 
M3= 0 
M5=M5+ 1 

IFIM5.LE.60) GO TO 150 
WR I TE < 6 » 20 10) 

M5= 1 

150 GO TO L, (160,170, 180) 


10 



160 ASSIGN 180 TO L 

WR I TE (6,2500) N6» < ( W ( J , K ) , J = 1 , 3 ) , K= 1 , N4 ) 

WRITE! 7,25 00) N6 , ( < W < J , K ) , J = 1 , 3 > , K= 1 , N4 ) 

M7 = l 

GO TO 190 

170 ASSIGN 180 TO L 

WRITE! 6,24 00) N6,N 6, ( ( W < J , K ) , J= 1 ,3 ) , K= 1 , N4 ) 

WRITE (7 ,2400) N6,N6, ( ( W ( J , K ) , J= 1 , 3 ) , K=1 , N4 ) 

M7 = 1 

GO TO 190 

180 WRITE! 6, 2 3 00) ( (W(J,K) ,J=1,3),K = 1,N4) 

WRI TE (7 ,2300) ( ( W ( J , K ) , J= 1 , 3 ) , K= 1 , N4) 

M7 =M7+ 1 

IF! M7.E0. 20) ASSIGN 170 TO L 
190 CONTINUE 
200 CONTINUE 

WRITE! 7,27 00) 

GO TO 100 
1000 FORMAT! I 1 ) 

2000 FORMAT ( 1H1»20( /) * 59X » ' LI TERAL FORM OF'// 

*52 X » ' T HE INVARIANT COEFFICIENTS GF • / / 

* 49 X , • THE CHARACTERISTIC POLYNOMIAL FOR A'/// 

* 64X , I 1, • X ', I 1/// 

* 64X, 'MATRIX ') 

2010 FORMAT! 1H1 ) 

2100 FORMAT (6X,'F(',I1,')=',4X, , +1') 

22 00 FORMAT! 6X, , F(',11,') = ',4X,8<A1, , A( ',11,',', 11,')')) 

2 300 FORMAT! 5X, •*',9X,8(A1,'A(',U, , ,',I1,')')) 

2400 FORMAT! 6X, , F(',I1,') = , ,'F(',11, , )',8(A1,'A(',I1,» ,' , I 1, • )')) 

2 500 FORMAT <6X,'F(',I1,') = , ,4X,8<A1, , A(',I1,',',I1,')')> 

999 STOP 

2600 FORMAT ( 'C '/ 'C' ,5X ,' COEFFICI ENT POLYNOMIALS FOR MATRIX OF ORDER', 

* 1 2 / ' C ' ) 

2700 FORMAT! 'C'/'C' /'C' ) 

END 
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SUBROUTINE COMB I N ( L t N t INDX ) 


PURPOSE 

TO PRODUCE THE NEXT COMBINATION (IN ASCENDING ORDER) OF L THINGS 
TAKEN N AT A TIME. NUMBER OF CALLS REQUIRED TO GET ALL 
COMBINATIONS IS BINOMIAL COEFFICIENT (L N) 


CALLING SEQUENCE C AL L C 0 MBI N ( L »N, I NDX ) WHERE 

L = THE NUMBER OF THINGS 

N = THE NUMBER OF THESE THINGS BEING COMBINED 

INDX = THE OUTPUT INDEX 

INDX ( 1 ) MUST BE ZERO ON THE FIRST CALL FOR A GIVEN N, L 
I NDX ( I ) WILL BE SET TO ZERO AFTER THE LAST COMBINATION 


OUTPUT 

THE OUTPUT IS SENT TO THE CALLING ROUTINE IN INDX. THE FIRST 
CALL WILL PRODUCE IN I NDX ( I ) TO I NDX ( N ) THE VALUES 1,2,...,N. 
SUCCESSIVE CALLS WILL PRODUCE, IN ASCENDING ORDER, THE OTHER 
COMBINATIONS OF THESE NUMBERS. BY ASCENDING ORDER IS MEANT, 
E.G., 1,2, 3, 4 IS LESS THAN 1,2, 3, 5 IS LESS THAN 1,2, A, 5 ETC. 


SAMPLE OUT PUT FOR 
AND 

CALL NO. 1 2 

INDX(l) 1 1 

I NDX ( 2 ) 2 2 

INDX ( 3 ) 3 A 


11 SUCCESSIVE CALLS TO COMBIN WITH L=6, N=3, 
I NDX ( 1 ) =0 ON THE FIRST CALL 
3 A 5 6 7 8 91011 

1 1 1 1 2 2 2 3 0 

233A33AAA 
5A55A5555 


RESTRICTIONS 

THE VALUES OF N AND L AND THE CONTENTS OF INDX SHOULD 
NOT BE CHANGED DURING THE (N L> CALLS TO COMBIN UNLESS A NEW 
SEQUENCE OF COMBINATIONS IS TO BE STARTED. IN THIS CASE, INDX(l) 
MUST BE RESET TO 0. 


DIMENSION I NDX ( 1) 

I F( INDXI 1 ).GT.O)GO TO 110 
DO 100 1=1, N 
100 I NDX ( I )=I 
RETURN 

110 I NDX (N)=INDX(N)+1 

I F( INDX ( N ) . LE .L ) RETURN 

L2=L 

M=N 

120 M=M-1 

I F( M.EQ.O )G0 TO 1A0 
L2 =L 2- 1 

INDX(M)=INDX(M)+1 
IF(INDX(M).GT.L2)G0 TO 120 
J=M+ 1 

DO 130 I = J ,N 
130 INDX ( I ) =INDX( I-D+l 
RETURN 

1 AO I NDX ( 1 ) =0 
RETURN 
END 
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SUBROUTINE PERMUT 
C PURPOSE 

C TO PERMUTE N ITEMS ON N FACTORIAL CALLS TO SUBROUTINE 

C 

C 

C METHOD 

C ON FIRST CALL TO PERMUT, IFIRST MUST EQUAL ONE, THE N ITEMS TO 

C BE PERMUTED MUST BE IN POSITIVE, ASCENDING ORDER IN 

C LAI 1,1) TO L A ( N, 1 ) , AND N MUST EQUAL THE NUMBER OF ITEMS BEING 

C PERMUTED. THIS FIRST CALL IS USED TO INITIALIZE COUNTERS AND 

C THE REST OF THE LA ARRAY, AND TO SET UP THE SIGN ARRAY (IS) WITH 

C POSITIVE VALUES. 

C 

C 

C OUTPUT 

C THE PERMUTED ITEMS WILL BE RETURNED TO THE CALLING PROGRAM IN 

C LAI 1,1) TO LA(N,1) WHILE THE SIGN OF THE PERMUTATION IS 

C RETURNED IN I SI 1 ) . THE SIGN IS INDICATED BY +1 OR -1 AS THE 

C PERMUTATION IS EVEN OR ODD. THE FIRST CALL TO PERMUT RETURNS 

C LAI 1,1) TO L A I N , 1 ) IN THE SAME ORDER AS IT WAS RECEIVED WHILE 

C THE SIGN IS RETURNED AS +1. HENCE THE ORIGINAL ORDER MUST BE 

C POSITIVE 

C 
C 

C SAMPLE OUTPUT FOR N= A , LA(1,10) TO LA(4,10) = 2,4,5, 7 
C 

C CALL NUMBER 1 23 4 5 67 8 9 10 11 12 

C LAI 1 ,1 ) 222222777777 

C LAI 2 , 1 ) 447 75 5225544 

C LAI 3,1) 57457 4452452 

C L A 1 4 , 1 ) 755447544225 

C ISI 1 ) +1 “I +1 “I + 1 -i "I + 1 “1 +1 

C 

C CALL NUMBER 13 14 15 16 17 18 19 20 21 22 23 24 

C LAI 1, 1) 555555444444 

C LAI 2,1) 774422552277 

C L A I 3 , 1 ) 247 247725725 

C LA! 4,1) 422774277552 

C I S I 1 ) +1 -1 +1 "I +1 -i -i + i -i "I +1 

C 

C 

C RESTRICTIONS 

C THE LA ARRAY, THE IS ARRAY, AND N MUST NOT BE CHANGED DURING THE 

C N FACTORIAL SUCCESIVE CALLS TO PERMUT UNLESS A NEW PERMUTATION 

C IS TO BE STARTED BY SETTING IFIRST BACK TO 1. IFIRST MUST NOT RE 

C CHANGED UNLESS A NEW PERMUTATION SET IS TO BE STARTED. 

C THE ITEMS TO BE PERMU TE D MUST ALL BE DIFFERENT AND, INITIALLY, 

C IN POSITIVE, ASCENDING ORDER. OTHERWISE, THE SIGN WILL BE WRONG. 

C 

C 

COMMON/NP/ I S I 9 ) ,LA( 9,9) ,N, IFIRST 
DIMENSION L I 9 ) , L L I 9 ) 

GO TO I 100 ,200 ) , IFIRS T 
100 I F I R ST = 2 
M = N-1 

DO 110 J =1 , M 
L I J ) = 0 
LLI J >=J 
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no i s ( j ) =1 

I F< N.EO.2 (RETURN 
DO 120 J=2,M 

DO 120 I =1 * N 
120 L A ( I , J)=LA (1,1) 

RETURN 

200 L ( 1 ) =L ( 1 > + 1 

IF( L (1) . GT. 1 )G0 TO 300 
K -L A (N * 1 ) 

LA( N, 1 ) = LA (N-l , 1) 

L A ( N- 1 , 1 ) =K 
IS( 1 ) =- 1 S ( 1 ) 

RETURN 
300 M = 1 
410 L ( M ) = 0 
M=M+ 1 

L ( M ) =L ( M ) + 1 

IF( L (M) .GT.LL(M) )G0 TO 410 
K =L A ( N ,M ) 

DO 420 1=1, M 

420 L A ( N- I + 1 ,M ) =LA ( N—I ,M) 

M2 =N— M 

L A ( M2 , M) = K 

M3=M-1 

I F(M0D(M,2 l.EO.l) I S(M)=-IS(M) 
DO 430 J =1 , M3 
I S ( J)=IS(M) 

DO 430 I=M2,N 
430 L A ( I , J ) = LA ( I , M ) 

RETURN 

END 
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